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In this talk, we present the arguments, that a new QCD regime - gluon saturation, has 
been reached at HERA. 



I. MAIN QUESTIONS AND PROBLEMS 

We hope that everybody agrees, we have two principle problems in low x physics: (i) matching between 
"soft" (npQCD) and "hard" (pQCD) processes; and (ii) theoretical description of high parton density QCD 
(hdQCD). My credo is [fj]-[|: these two problems are correlated and the system of partons always passes the 
stage of hdQCD ( at shorter distances ) before it goes to the black box, which we call non-perturbative QCD 
and which, in practice, we describe in old fashion Reggeon phenomenology. 
In this talk you can find the answer to the following questions: 

• Is there any violation of the DGLAP evolution? 

• Can we describe matching between "soft" and "hard" processes? 

• Do we see a signal of the gluon saturation in HERA data? 

• What is a current situation in theory for high parton density QCD system? 

Unfortunately, the lack of space does not allow us to discuss such a hot problem as the status of the BFKL 
Pomeron as well as the manifestation of high parton density QCD in the Tevatron data. 

II. TWO SCALES OF DIS 

About thirty years ago Gribov Q noticed that photon - hadron interaction at high energies has two distinct 
stages as far as time - space picture of interaction is concerned: 
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FIG. 1. 



1. 7* — ► hadron system ( qq - pair ); 



2. hadron system ( qq - pair ) interacts with the target. 
Therefore, we can describe photon-hadron interaction as follows 

Statin* P) = l*»l 2 a tot(np) , (1) 

n 

where n denotes the set of quantum numbers which diagonalize the high energy interaction matrix. This 
set of quantum numbers we call the correct degrees of freedom (DOF) and <Jtot{np) is the total cross section 
for the interaction of the hadron or parton system with quantum numbers n with the target. Eq. (|l|) is 
useful only if we know the correct DOF. Fortunately, we do know them at short distances where n are colour 
dipoles U and Eq. (Q) reads as 

vtattfp) = J d 2 r t £ dz\*(Q 2 ;r t ,z)\ 2 a tot (rlx) . (2) 



A. Separation scale r^ p = 1/M 



However, at long distances we do not know the correct DOF. The scale which says what distances are 
short we call a separation scale (see Table. 1). 
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Table 1. 

, sen , 

r± < r ± y < r± 

^(Q 2 ;r tj z)^ pQCD ^(Q 2 ;r t ,z)^ npQCD 

DOF: colour dipoles || DOF: constituent quarks | 

colour dipoles (7j 

cr to t (r± , x) oc cr £o£ (r_L , VP) -> 

xG{x, 4/r^_) Regge phenomenology 



B. Saturation scale r^_ at « 1/Q s (x) 

This scale can be estimated from the equation 

_ 3 7r 2 a s xG(x,Q 2 s (x)) _ 
K ~ 2Ql(x) ttR 2 ~ lj (3) 

which says that the packing factor of partons in the parton cascade is about unity or, in other words, at this 
scale the parton system is so dense that we cannot apply the standard methods of perturbative QCD. 

The physical meaning of ft is clear from Fig. 2 which shows the parton distributions in the transverse 
plane. At short distances ft is the product of the parton density in transverse plane which is equal to 
number of partons /area = xG(x, Q 2 )/ttR 2 multiplied by the cross section of the parton interaction oc as/Q 2 . 
Therefore, at the scale where ft « 1 the interaction between partons become essential and new approach 
should be developed in QCD to describe the high parton density system. 

The natural hierarchy of the scales is r^_ a£ <C r s ^ p accordingly to our main idea. 



III. THEORY STATUS 



The theory of high density system in QCD is in a very good shape now. The problem has been attacked 
from two different point of view: (i) from the pQCD region by summing corrections to the evolution 

equations due to high density of partons; and (ii) from the non-perturbative QCD region by developing 
effective Lagrangian approach [||J^] dealing with such a system. The resulting evolution equation that has 
been proven [|Io| looks as follows 

a el (xoub t ,y) + — / 

77 J p X 02 X 12 

,2/) - a el (x 2:b u y) a el (x 12 ,b t ,y) ) , (4) 

where a el is the elastic amplitude for dipole scattering at fixed impact parameter b t and at energy s ( 
y = Ins). 



da el (xqi, b u y) _ _ 2C F a s ( x\ x \ 
dy 7T V P 2 J 

■ ( 2a el (x 02 ,bt - ^x 12 
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FIG. 3. 

The physical meaning of Eq. (0) is very clear from Fig. 3. Indeed, with probability 2 Xoi 2 dipole with 

M_ 1 X 02 X 12 

size Xq X decays in two dipoles with sizes Xq 2 and x^ 2 (see Fig. 3 ). These two dipoles interact with the 
target: each produced dipole interacts with the target separately ( the second term in Eq. (0) ) or two 
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dipoles interact with the target simultaneously ( the third term in Eq. (|4|) ). The first term describes the 
fact that dipole Xq X disappears from the initial state after decay into two dipoles. This equation, which 
has been suggested in the momentum representation in Ref. [Q, has a lot of nice properties including the 
correct matching with the DGLAP evolution equations. However, the most important message concerning 
this equation is that this equation can be derived using the effective Lagrangian approach j|. It should be 
stressed that we know not only the equation but also the initial condition for it. 



IV. MATCHING OF "SOFT" AND "HARD" PHOTON-PROTON INTERACTIONS 



Serious attempts {H) to find the value of the separation scale r s _^ p = 1/Mo were undertaken using 
Gribov's formula || during the past decade, starting from pioneering paper of Kwiecinski and Badelek fl^] . 
Gribov's formula reads (see Fig. 5) 



17 J 3^ J (Q 2 + M 2 ) 1 ' ' } (Q 2 



(Q 2 + M' 2 ) 




FIG. 4. 

The following assumptions are made to describe the matching between "soft" and "hard" processes: 

M < M) M > M 

M = M' M ^ M' 

D-L Pomeron + AQM "Hard" Pomeron = pQCD approach 

The result ( see Fig. 5 ) is that the separation scale depends on polarization of the incoming photon and 
it is equal 0.7 < Mq T < 0.9 GeV 2 for transverse polarized photon and it is smaller for the longitudinal 
polarized photon ( M$ L < 0.4 GeV 2 ). 



V. WHERE ARE SC? 



The HERA data put a puzzling question on the table: 

• On one hand, the data can be described by the routine DGLAP evolution equations 

• On the other hand, the gluon density measured at HERA is so high that some effect of hdQCD 
should be seen (see Fig. 6 where parameter n is plotted as it appears in HERA data. ) 
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FIG. 5. An example of description the HERA experimental data using Gribov's formula, taken from GLMN paper 



The natural question to ask is where to look for a saturation scale Q s (x) . Fig. 7 shows us that the 
high parton density corrections to F2 is rather small while they are substantial for the gluon structure 
function. 



VI. Q|(X) FROM Q 2 - DEPENDENCE OF F 2 - SLOPE 



There is a hope that the i^-slope ( q\^q2 ) will provide a measurement of the saturation scale since in 
the DGLAP evolution this slope 9 ^q 2 = i^f x G DGLAP (x, Q 2 ) directly proportional to the gluon 
structure function. On the other hand, a gluon saturation leads to the slope which is proportional [Q 
to Q 2 R 2 at fixed x where R is the target size. Indeed, it turns out that hdQCD corrections are 
able to describe all experimental data on the F2-slope (see Fig. 8). We consider as an important sign, 
that we are on the right track, the fact that two DGLAP parameterizations GRV94 and GRV'98 lead 
to a good description of the experimental data after taking into account hdQCD effects. 
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FIG. 6. 
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[Results:] 

► SC are large for 
xG(x,Q 2 ) but their 
values do not depend on 
the way how we take SC 
into account; 

► SC are rather small for 

^ 2 ; 





FIG. 7. The hdQCD corrections to Fi and xG(x,Q 2 ) calculated in Ref. |l5|. Calculation were done using Eq. (| 




FIG. 8. 

However, the i^-slope data can be equally well described by "soft" and "hard" contribution without 
hdQCD effects if we assume that the "soft" contribution stems from rather short distance 0.3 0.5 fm 
(see Fig. 9). 



VII. HIGH DENSITY QCD AND DIFFRACTIVE J/^-PRODUCTION IN DIS 

It turns out that it is very instructive to consider two observables : the F2- slope and the J/^-production 
in DIS. The hdQCD effect are essential in both observables and a simultaneous analysis of them could give 
an information on the value of the saturation scale. Performing such an analysis we obtain the following: 
(i)none of MRS parameterizations survives; (ii)all GRV parameterizations survive only with hdQCD effects; 
and (ii) x 2 /n.d.f. is excellent for GRV + hdQCD effects (see Fig. 10 ). 
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FIG. 9. The F2 -slope in D-L model with "soft" contribution at rather short distances 
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FIG. 10. J/\Ij production in hdQCD approach. Curves are taken from Ref. jnj 



VIII. A NEW SCALING 



The saturation hypothesis got a new support by Golec-Biernat and WiisthofT [[L8| who suggested a simple 
model that introduces only one scale for colour dipole - target interaction: the saturation scale Q s {x) , 
which they defined by fitting experimental data. This simple model does not take into account even such 
well established property as scaling violation but describes all experimental data from HERA for x < 0.01. 

In Ref. [H9 it was noticed that the HERA data show a new scaling: cr(j*p) = R 2 F(Q 2 /Q 2 (x)) for all 
values of Q 1 ^ at x < 0.01 ( see Fig. 11 ). Such a scaling was expected |i|,^o|,p|,^H in the saturation region 
to the left of the critical line k = 1 in Fig. 2. The value of the saturation scale Q s {x) can be measured as 
a value of Q 2 at which we see a deviation from this scaling. The fact, that all data at HERA kinematic 
region for x < 0.01 show this scaling, confirms the idea that at x < 0.01 DIS is deeply in the saturation 
region. It should be stressed that the whole idea of saturation is the simple fact that the only one scale 
f saturation determines the scattering amplitude for the distances longer than r satur ation' This very fact one 
can see directly in Fig. 2, noticing that a hadron looks as the diffraction grid with the size v 'saturation ifr 

the 

saturation region (ft ^> 1 in Fig. 2 ). 




a(y*p) = F(T=Q/Q(x)) 

• Bartles & E.L. (1992) ; McLerran & Venugopalan (94) 
E.L. & Tuchin (1999) ;McLerran & Kovchegov (98); 

FIG. 11. 
IX. SUMMARY 

The brief review on low x physics, which is given in this talk, allows us to make a definite conclusion, 
that a new QCD regime has been reached at HERA: the regime of high parton density QCD with a gluon 
saturation. Let us list all arguments for such a new regime: 

• There are no experimental data that are in a contradiction with the asymptotic prediction of high 
density QCD. Indeed, 

- The F 2 -slope shows dF 2 /dlnQ 2 oc Q 2 behaviour at Q 2 < Q 2 (x) |§|,|j§; 

— The ratio of diffraction cross section in DIS to the total DIS cross section is constant at HERA 
kinematic region |2l|| ; 
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The inclusive diffraction stems from short distances as should be in a gluon saturation picture 
in which a hadron looks as a diffraction grid with a typical size 1/Q 8 (x) [p4|,p5|; 



• The HERA data show a new scaling |19|, predicted theoretically [|l|,|20|,pU21f| in the saturation region 
, namely, that a(>y*p) = R 2 F (Q 2 / Q 2 s (x)) for all values of Q 2 at x < 0.01. 

• All data can be described in the simplest saturation model of Golec-Biernat and Wiisthoff fll8|j ; 

• The current parameterizations, based on the DGLAP evolution equations, cannot describe simul- 
taneously the Q 2 -behaviour of the F2 slope and energy behaviour of the J/\I/ production in DIS and 
photoproduction |T^ ]; 

• The theoretical developments in high density QCD has been so remarkable during the past two 
years, that we can trust Eq. (|4|) which predicts the essential high density collective phenomena in 
HERA kinematic region; 

• The simple estimates based on the gluon density extracted from the HERA data show large packing 
factor (see Fig. 6 ) or, in other words, they show the strong shadowing corrections are needed in HERA 
kinematic region. 

In spite of everything mentioned above, not everybody will agree with the strong statement that has been 
formulated here. The reason for this is very simple: most of data, that has been considered above, have 
different explanations without SC. All these alternative explanations cannot be considered as natural ones, 
but it is behind many people scepticism on a new QCD regime at HERA, that the DGLAP evolution has 
more fundamental origin in QCD than all estimates with SC. It is not true at all and SC are more general 
approach than the DGLAP equation because (i) they are consistent with the s-channel unitarity; and (ii) 
the non-linear equation (see Eq. (|J) )has the same deep operator proof as the DGLAP evolution equation. 

Acknowledgments : 

I wish to thank E. Gotsman, M. Lublinsky, U. Maor, E. Naftali and K. Tuchin for a pleasure to work with 
them and to discuss with them all problems of low x physics. 

This research was supported by part by Israel Academy of Science and Humanities and by BSF grant # 
98000276. 



[1] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep 100 (1983) 1. 
[2] A.H. Mueller and J. Qiu, Nucl. Phys. B268 (1986) 427. 

[3] L. McLerran and R. Venugopalan,P%s. Rev. D49 (1994) 2233,3352, 50 (1994) 2225, 53 (1996) 458, 59 (1999) 
094002. 

[4] V.N. Gribov, Sov. Phys. JETP 30 (1970) 709. 
[5] A.H. Mueller, Nucl. Phys. B415 (1994) 373. 

[6] E. Levin and L. Frankfurt, JETP Lett. 2 (1965) 65; H. J. Lipkin and F. Scheck, Phys. Rev. Lett. 16 (1966) 71. 
[7] J. R. Forshaw, G. Kerley and G. Shaw, Phys. Rev. D60 (1999) 074012, Nucl. Phys. A675 (2000) 80; 

M. McDermott, L. Frankfurt, V. Guzey and M. Strikman, "Unitarity and the QCD-improved dipole picture", 

hep-ph/9912547 . 
[8] E. Levin and M.G. Ryskin, Phys. Rep. 189 (19267) 1990; 

J.C.Collins and J. Kwiecinski, Nucl. Phys. B335 (1990) 89; 

J. Bartels, J. Blumlein and G. Shuler, Z. Phys. C50 (1991) 91; 

E. Laenen and E. Levin, Ann. Rev. Nucl. Part. Sci. 44 (1994) 199 and references therein; 
A.L. Ayala, M.B. Gay Ducati and E.M. Levin, Nucl. Phys. B493 (1997) 305, B510 (1990) 355; 
la. Balitsky, Nucl.Phys. B463 (1996) 99; 

Yu. Kovchegov, Phys. Rev. D54 (191996) 5463, D55(1997) 5445, D60(2000) 034008, D61 (2000)074018; 
A.H. Mueller, Nucl. Phys. B572 (2000) 227, B558 (1999) 285; 
Yu. V. Kovchegov, A.H. Mueller, Nucl. Phys. B529 (1998) 451; 
E. Levin and K. Tuchin, Nucl. Phys. B573(2000) 833. 



11 



[9] J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert, Phys. Rev. DD55 (1997) 5414; 
J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. D59 (1999) 014015; 
J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. D59 (1999) 014015; 

J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Phys. Rev. D59 (1999) 014014,034007, Erratum- 
ibid. Phys. Rev. D59 (1999) 099903; 

A. Kovner, J.Guilherme Milhano an d H. Weigert, OU TP-00-10P 1 NORDITA-2000-14-HE 1 |hep-ph/0004014 ; 

H. Weigert, NORDITA-2000-34-HE, |hep-ph/0004044 
[10] la. Balitsky, Nucl.Phys. B463 (1996) 99; Yu. Kovchegov, Phys. Rev. D60 (191000) 034008. 
[11] E. Gotsman, E. Levin and U. Maor, Eur. Phys. J. C5 (1998) 303; A.D. Martin, M.G. Ryskin and A.M. Stasto, 

Eur. Phys. J. C7 (1999) 643; E. Gotsman, E. Levin, U. Maor and E. Naftali, Eur. Phys. J. CIO (1999) 689. 
[12] B. Badelek and J. Kwiecinski, Phys. Lett. B295 (1992) 263. 

[13] A. M. Cooper-Sarkar, R. C. E. Devenish and A. De Roeck, Int. J. Mod. Phys. A13 (1998) 33; 

H. Abramowicz and A. Caldwell, Rev. Mod. Phys. 71 (1999) 1275. 
[14] E. Gotsman, E. Levin, U. Maor and E. Naftali, Nucl.Phys. B539 (1999) 535; E. Gotsman, E. Levin and U. 

Maor, Phys. Lett. B425 (1998) 369. 
[15] A.L. Ayala, M.B. Gay Ducati and E.M. Levin, Nucl. Phys. B493 (1997) 305, B510 (1990) 355. 
[16] A. Donnachie and P.V. LandshofT: Phys. Lett. B437 (1998) 408; Phys. Lett. B470 (1999) 243. 
[17] E. Gotsman, E. Ferreira, E. Levin, U. Maor and E. Naftali, " Screening corrections in DIS at low Q 2 and x" , 
Talk given at 30th Internat ional Conference on High-Energy Physics (ICHEP 2000), Osaka, Japan, 27 Jul - 2 
Aug 2000, |hep-ph/0007274[ 



[18] K. Golec-Biernat and M. Wusthoff, Phys. Rev. D59 (1999) 014017; D60 (1999) 114023; K.Golec-Biernat,Talk 
at 8th In ternational Work shop on Deep Inelastic Scattering and QCD (DIS 2000), Liverpool, England, 25-30 
Apr 2000 jiep-ph/000608C| . 

[19] A.M. Stas to, K. Golec-Bier nat and J. Kwiecinski, " Geometric scaling for the total j*p cross-section in the low 



x region", [hep-ph/0007192 



[20] J. Bartels and E. Levin, Nucl. Phys. B387 (1992) 617. 

[21] Yu.V. Kovchegov and L. McLerran, Phys. Rev. D60 (1999) 054025. 

[22] B. Foster (ZEUS Collaboration), Invited talk, Royal Society Meeting, London, May 2000; 

ZEUS collaboration: J. Breitweg et al., Phys. Lett. B487 (2000) 53. 
[23] M. Klein (HI collaboration), "Structure Fun ctions in Deep Ine lastic Lepton-Nucleon Scattering" Talk at Lepton- 

Photon Symposium, Stanford, August 1999, |hep-ex/0001059| . 
[24] ZEUS Collaboration: J. Breitweq et al., Eur. Phys. J. C6 (1999) 43. 
[25] HI Collaboration: C. Adloff et al., Z. Phys. C76 (1997) 613. 



12 



